ON NEWTON INTERPOLATION OF SYMMETRIC 

FUNCTIONS. A CHARACTERIZATION OF 
INTERPOLATION MACDONALD POLYNOMIALS. 



0\ ' Andrei Okounkov 



(^ 



Contents 



Q ■ 1 . Brief introduction 

«^ . 2. General interpolation problem 

CN I 3. Examples of interpolation polynomials 
3.1 Universal interpolation polynomials 

>- ■ 3.2 Factorial monomial symmetric functions 

^ ■ 3.3 Factorial Schur functions 

Q ! 3.4 Interpolation Macdonald polynomials 

CN [ 4. Statement of the characterization theorem 

5. Reductions of the proof 

0\ • 6. Proof of the theorem in the two variables case {n = 1) 

bJO' 7. Proof of the theorem for n > 1 

o3 ! A. Appendix. Table of perfect grids 

qh| B. Appendix. An example of the universal interpolation polyno- 

IJ ' mial 

X 

1. Brief introduction 

The purpose of this paper is to characterize interpolation Macdonald polynomials 
inside a very general Newton interpolation scheme for symmetric polynomials. This 
general Newton interpolation problem is discussed in Section 2; it depends as on a 
parameter on a map fl 

^Z^ofoi--ri.hies_J]_^ ground field k, 

which we call a grid in k. Our present understanding of this general problem can be 
described as follows: most of it is covered by an unexplored and mysterious ocean 
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2 ANDREI OKOUNKOV 

formed by generic grids O (see Section 3.1). In the midst of this abyss there are 3 
pieces of dry land, namely the 3 following exactly solvable cases: 

(1) factorial monomial symmetric functions, 

(2) factorial Schur functions, 

(3) interpolation Macdonald polynomials, 

which are described in Sections 3.2-3.4. The first case is just dull, the second one is 
still rather elementary; both of them are parameterized by an arbitrary sequence of 
pairwise distinct elements of the ground field k. These two continents are joined by 
a beautiful archipelago of interpolation Macdonald polynomials. More precisely, by 
interpolation Macdonald polynomials we mean the so called SC-type interpolation 
Macdonald polynomials, introduced and studied in [Ok4]. As particular cases and 
degenerations these polynomials include polynomials studied by F. Knop, G. 01- 
shanski, S. Sahi, and the author in a long series of papers, see References. These 
polynomials depends on 5 parameters of which only 3 are non-trivial because of an 
action of a 2-dimensional group of affine transformations. Some of the properties 
of these most remarkable polynomials are discussed in Section 3.4. 

It is natural to ask if any simple abstract property characterizes the 3 above 
exactly solvable cases of our general interpolation problem. As such a property we 
propose the extra vanishing property (4.2) which says that the Newton interpolation 
polynomials should vanish not only at those points where they are supposed to 
vanish by their definition but also at certain extra points "for free" . More precisely, 
the polynomial labeled by a partition ^ vanishes at the point labeled by a partition 
A unless // is less or equal to A in the partial order of partitions by inclusion // C A . 
We call all grids that enjoy this property perfect. 

This extra vanishing property can be compared to the following well-known 
property of ordinary Macdonald polynomials. Although the Gram-Schmidt or- 
thogonalization process requires a choice of a total order on the polynomials to be 
orthogonalized, the Macdonald orthogonal polynomials do not actually depend on 
the choice of a total order on the monomial symmetric functions as long this total 
order is compatible with the partial dominant order of partitions.^ 

Also, the extra vanishing property can be compared to a well known phenomenon 
in integrable systems where many exactly solvable systems have "extra" integrals 
of motion, that is, more integrals of motion than is required by the definition of 
integrability [Kr,CV]. It is interesting to notice that certain integrable many-body 
systems to which interpolation Macdonald polynomials are very closely connected 
were conjectured to be characterized by this "extra" integrability property [CV]; 
later, however, certain new examples were found in [CFV]. 

Our situation is, of course, much simpler. Our main result (Section 4) is that the 
three above cases plus degenerations of the third one exhaust the set of all perfect 
grids. As a corollary of this theorem, we conclude that no other grid O admits a 



^It would be probably interesting to describe all interpolation or orthogonal symmetric poly- 
nomials which satisfy such a "extra triangularity" condition. 
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tableaiix sum formula of the form (4.1) for the Newton interpolation polynomials, 
nor does it admit an integral representation of interpolation polynomials analogous 
to the g-integral representation for interpolation Macdonald polynomials obtained 
in [Ok4] . That is, any new exactly solvable case of symmetric Newton interpolation 
has to be based on some entirely new type of formulas. 

The proof of this characterization theorem is given in Sections 5-7. Section 5 
contains some general statements, whereas the two other sections are devoted to 
the consideration of the many possible cases. 

There exist also non-symmetric Macdonald interpolation polynomials which form 
a linear basis in the algebra of all polynomials, see [Kn,S2]. It is plausible that 
those polynomials might have a similar characterization. Note also that a certain 
characterization of ordinary Macdonald polynomials inside some general class of 
orthogonal polynomials was found by S. Kerov in [K]. 

2. General interpolation problem 

We consider Newton interpolation of symmetric polynomials in n + 1 variables 
(n = 0, 1, . . . , oo) with coefficients in some infinite field k. First assume for simplic- 
ity that n < oo; the case n = oo will be covered at the end of this section. Any 
natural basis in the space of symmetric polynomials of degree < d is indexed by 
partitions^ 

A = (Ao > Ai > ■ ■ ■ > An > 0) 

of length £{X) < n + 1 (recall that i{X) is the number of non-zero parts in A) such 
that 

|A| <d. 

Here |A| = Aq + ■ ■ ■ + A„. Therefore, the knots of our Newton interpolation should 
be also indexed by partitions. In other words, we need a function 



O: {A,£(A) <n + l} 



,n+l 



which takes a partition A to the corresponding knot of interpolation 0(A). We can 
construct such a function in the following way. Choose a function 

(2.1) O : {0, ...,?i} X Z>o — >k, 

which we shall call a grid in k, and then set 

0(A) := (0(0, Ao), . . . , 0(n, A^)) G k'^+i . 

To simplify notation, we shall use abbreviations 

[^j]:=0(z,j), A:=0(A). 



■^Our notation conventions about the partitions are slightly different from the standard ones 
used in [Ml]. For example, we number the parts of partitions as well as coordinates of any vector 
starting from zero. This leads also to a different definition of a diagram of a partition but makes 
certain formulas look more symmetric. 
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Definition 2.1. A grid O is said to be non- degenerate if for any partition /x of 
length < n + 1 there exists a symmetric polynomial 



P^(xo, . . . , Xn; O) G k[xo, . . . , X 



7 "^nj 



lS(n+l) 



satisfying the following Newton interpolation conditions 

(1) the degree of Pfj.{x; O) is < \fx\; 

(2) -P^(A; O) = for all partitions A of length < n + 1 such that |A| < |//| and 

(3) P^(/2;O)^0. 

Remark 2.2. It is clear that if the grid O is non-degenerate then 

(1) all polynomials P^{x]Q) are uniquely defined up to a non-zero factor from 
the field k; 

(2) the polynomials P^{x; O) as fu, ranges over all partitions of length < n + 1 
form a linear basis of the vector space k[xo, . . . , Xn]^^^~^^' of all symmetric 
polynomials; 

(3) the degree of P^{x; O) is precisely \n\. 

The Newton interpolation polynomials P^{x]Q) generalize polynomials consid- 
ered by S. Sahi in [SI]. 

All non-degenerate grids admit a simple description. 

Proposition 2.3. A grid O is non- degenerate if and only if 

(2.2) [ij]^[i'j'], ^hj i>i',j<j'. 

Before we prove this proposition let us introduce the following operation on grids 
(2.1). For any m = 0, . . . , n we can define a grid 

flrn ■ {0, . . . , m} X Z>o — > k 

by simply restriction 0^(z, j) := 0(z, j) of the grid O. 

Now the following proposition follows immediately from the definition of poly- 
nomials P^{x; Vt). (Since so far these polynomials were defined up to a scalar factor 
only the equality (2.3) is to be understood for the moment as an equality up to a 
scalar factor. However, below in Definition 2.8 we shall choose a particular normal- 
ization of the polynomials P^{x; O) which will make (2.3) into a precise equality.) 

Proposition 2.4. Suppose Q is a non- degenerate grid and n is a partition such 
that fXn = 0- Then 

(2.3) P^{xo,...,Xn-i, [nO];0) = P(^o,...,^„_^)(xo, . . . , a;n-i; ^n-i) ■ 
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Corollary 2.5. // a grid O is non- degenerate then so are all grids Vtrn, ^ — 
0, ...,n. 

Proof of Proposition 2. 3. Suppose that Vt is non-degenerate and show that 

By the above Corollary we can assume that i = n. Suppose that, on the contrary, 
[nj] — [i'j'] and j' > j. Then the following polynomial 

n J -1 
m=0 fc=0 

is symmetric, has degree nj' and vanishes at all points A such that |A| < nj'. Thus, 
the "only if" part of the proposition is established. 

The "if" part of the proposition follows from the following argument, which is 
an abstract form of the argument of S. Sahi, see [SI]. 

Proposition 2.6. Suppose a grid O satisfies the conditions (2.2) and let S Ck be 
any ring such that 

m - Wj] 

Then for any d = 0,1, . . . and any function 

(2.4) (j): {X,\X\<d,i{X)<n + l} — >§ 
there exists unique symmetric polynomial 

f^ e S[xq,...,x 
of degree d such that 

(2.5) ^(A) = </>(A) . 

Proof. Induct on n and d. The cases n = or (i = are clear. Suppose n,d > Q. 
Introduce a linear map 

ext : §[a:o, . . . , x^^.^f^^^ -^ §[^0, . . . , xj^("+i) , 

which, by definition, takes a monomial symmetric function in variables 

xq - [n%...,Xn-i - [nO] 



S(n+1) 

, . . . , vt-TlJ 
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to the same monomial symmetric function in variables 

xo - [nO],...,Xn - [nO] . 
Observe that it is a degree preserving injection and that 

(2.6) {ext f){xo, . . . , Xn-i,[nO]) = /(xq, • • • , x„_i) . 
Now given (p we shall look for the solution /^ in the form 

n 

(2.7) /^ = ext/i + /2n(^^-H)' 
where /i and /2 are unknown polynomials such that 

/i G§[xo,...,x^_l]^('"^ deg/i=rf, 
/2G§[xo,...,a:,]^^"+^^ deg/2=rf-n. 

First, consider the equations (2.5) for partitions A such that An = 0. Since then 
the second summand in (2.7) vanishes these equations by (2.6) and inductive as- 
sumption determine the polynomial /i. 

Now consider the equations (2.5) for the remaining partitions A (that is, for A 
such that An > 0). Rewrite them in the form 

and observe that the RHS of (2.8) lies in §. The set of partitions A such that |A| < d 
and An > is in bijection 

Ah^A-i = (Ao-l,...,An-l) 

with the set of partitions such that |A| < d — n. Therefore, replacing O by O^ and 
using the inductive hypothesis we can find a polynomial /2 of degree d — n satisfying 
the equations (2.8). This proves existence of the required polynomial /<^. 

To verify uniqueness of fcj, it suffices to consider the case S = k. We have an 
obvious k-linear map from the space of symmetric polynomials / of degree < d to 
the space of functions (p of the form (2.4), namely 

f^<i>, <^(A):=/(A). 

Both spaces have same dimension and we just have proved that this map is surjec- 
tive. Hence it is an isomorphism with the inverse map 
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This concludes the proof of the Propositions 2.6 and 2.3. D 

Let us introduce another operation on grids (2.1). Given a grid O, we can for 
any /c = 0, 1, 2, . . . define a new grid O'^ 

O^ : {0, ...,n} X Z>o — >k 

by the formula 

0'=(z,j):=0(z,j + fc). 

Then the following property follows immediately from Definition 2.1 and Proposi- 
tion 2.3 (As in the case of Proposition 2.4, observe that since the normalization of 
the polynomials P^{x; Q) is yet to be specified the equality (2.9) is to be understood 
for the moment as an equality up to a scalar factor.) 

Proposition 2.7. Suppose Q is a non- degenerate grid. Then the grid Q^ is also 
non- degenerate and for any partition // such that Hn > we have 

n 

(2.9) P^(a:o,...,Xn;0) = P^_j(xo, . . . ,Xn;0^) W{xi - [nO]) , 

where /U — 1 stands for partition {fxo — 1, . . . , /Utj, — 1). 

Now suppose that O is a non-degenerate grid and let us specify the normalization 
of the interpolation polynomials P^{x; O). Let us identify any partition A with its 
diagram which is, by definition, the following subset of Z>o x Z>o 

A = {{i,3)^3 < Aj - 1} C Z>o X Z>o . 

Note that this definition differs by a coordinate shift by 1 from the standard defi- 
nition used in [Ml]. Let A' stand for the partition corresponding to the transposed 
diagram of A, that is 

(2.10) A;:=#{z,A. >j}. 

We are free to normalize the polynomials P^(x;0) by setting their value at the 
point /i to any non-zero element of k. We make the following choice 

Definition 2.8. Suppose O is a non-degenerate grid and /U is a partition with 
<n-\-l parts. Then we define P^{x; Q) to be the unique polynomial satisfying the 
conditions of Definition 2.1 and the following normalization condition: 

(2.11) P^(/2;0)= H ([i,^,]-[^;.-l,_^-]). 

(*J)eM 

The reason we prefer the normalization (2.11) is the following 
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Proposition 2.9. The normalization (2.11) is the unique normalization compatible 
with (2.3) and (2.9). 

Proof. Induct on n and |/u|. Since for any partition // either //^ > or //^i = we 
can apply either (2.3) or (2.9) to the evaluation of 

and thus reduce it to the case of smaller values of |/u| or n respectively. D 
By definition, let A C k be the subring generated by the elements 

1 



(2.12) A = Z 



M^JT. 



il Al\ 



[ij] - [i'j 
Then since by (2.11) we have 



Ck, i>i'j<j' 



for any fi we conclude from Proposition 2.6 and Definition 2.8 that 

Proposition 2.10. P^{x; O) G A[xo, . . . , Xnf^''+^^ . 

Remark 2.11. Now let us explain the definition of the polynomials P^{x; O) in the 
case of infinitely many variables (n = oo). Suppose we are given a grid 

(2.13) 0:Z>oxZ>o — >k. 

Then we need first a suitable definition of the algebra of symmetric polynomials in 
infinitely many variables. For any to = 0, 1, 2, . . . define a homomorphism 

(2.14) ResZ-^^ : k[xo, . . . , x^+i]^('"+2) -^ k[a:o, • • • , Xmf^'^^^'^ 
by the formula 

(Res^+V) ixQ,...,x^) = f{xQ,...,Xm, [m + 1,0]). 
Observe that 



m 



degRes:^+V<cieg/. 
By definition, let A" be 

A":=limk[a:o,...,a:J^(-+i) 



7 "^m.J 



the inverse limit of filtered (by the degree of polynomials) algebras with respect to 
homomorphisms (2.14). Observe that by construction for any 

/e A" 
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its degree deg / is well defined and so are its values at the points of the form 

/(A)=/([OAo],[lAi],[2A2],...), 

where A is a partition. Therefore the interpolation problem described in Definition 
2.1 is well defined for polynomials in the algebra A*^, so it makes sense to ask 
whether the grid O is non-degenerate and what are the corresponding Newton 
interpolation polynomials. Using definitions and Propositions 2.3 and 2.4 one easily 
checks that 

(1) a grid (2.13) is non-degenerate if and only if the conditions (2.2) are satisfied 
(which is also equivalent to non-degeneracy of all grids O^, < to < oo); 

(2) for any partition // the sequence 






defines an element of A^ which equals P^,{x] O). 

3. Examples of interpolation polynomials 

3.1 Universal interpolation polynomials. 

Let k" be the field 

k" = Q(w.,), (z,j)e{0,...,n}xZ>o, 
of rational functions in variables Uij and let the grid O" be the grid 

[iJV = Uij , 

which is clearly non-degenerate (the superscript u stands here for "universal" ) . We 
have the corresponding polynomials 

(3.1) P^{x;n^)eA^[xo....,x^f^^+'\ 

where the ring A" is, according to the definition (2.12), the ring 

1 



Z 



Uij 1 



Uij Ui'j' 



ck^ i>i'j<j'. 



The polynomials (3.1) are universal in the following sense. For any field k and any 
non-degenerate grid 

O : {0, . . . , n} X Z>o — > k 

we have a natural homomorphism of rings 
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which is well defined by non- degeneracy of O and takes universal interpolation 
polynomials to the corresponding specific ones 

One can compute a few universal polynomials explicitly (see an Example in 
Appendix B) but the formulas get very soon very complicated. One is forced 
therefore to look at least general but nicer examples. 

The three other examples we shall consider in this section will be exactly solvable 
in the following sense: there exist rather simple closed formulas for the interpolation 
polynomials. For simplicity we consider the case of finitely many variables (n < oo). 
However, all formulas are stable in the sense described in Remark 2.11. 

3.2 Factorial monomial symmetric functions. 

Suppose that 

[ij] = Cj , 

where cq, ci, . . . are pairwise distinct elements of k. Then it is easy to see that the 
interpolation polynomials P^ are simply factorial monomial symmetric functions 

^ n Hi — 1 

where ^^ stab^^^+i) n is the number of permutations s G S{n + 1) which leave the 
vector n = (fiQ, . . . , Hn) invariant. It is also easy to see that the polynomials (3.2) 
vanish at more more point then it is prescribed by their definition. Namely, we 
have 

P^(A;O)=0, unless /U C A , 

where the notation // C A means that Hi < Xi for all i (which is equivalent to the 
diagram of fx being a subset of the diagram of A). 

3.3 Factorial Schur functions. 

Suppose that 

M = Cj-^ , 

where . . . , c_i, cq, ci, . . . are pairwise distinct elements of k. Introduce the following 
factorial Schur polynomial (see [M2] and references therein and also [001]) 



det \{Xi — C-n) ■ ■ ■ (Xi — Cu-_o_i I I „ ,. ./ 

(3.3) s,(x; Q) := ^' "^ '. "'' , '^^'"■-^" • 

ll0<i<j<nl"^* "^JV 

First observe that (3.3) is indeed a polynomial because the numerator in (3.3) is an 
anti-symmetric polynomial in xq, . . . , a:^ and hence is divisible by the denominator. 
It is also clear that the ratio has degree \\x\ and its top-degree component is the 
classical Schur function s^(x). 
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We shall check momentarily (borrowing the argument from [Okl], Section 2.4; 
see also [001]) that we have 

(3.4) P^(x;0) = s^(a;;0). 
First show that 

(3.5) s^(A;O) = 0, unless /U C A . 

Observe that the denominator in (3.3) does not vanish at any point of the form A. 
Therefore, it suffices to check that the numerator vanishes unless // C A. Suppose 
that Afc < Hk for some k. Then for all < j < /c < z < n we have 

and the corresponding matrix element vanishes 

(Cxi-i - C-n) ■ ■ ■ {Cxi-i - C^.-j-i) = . 

Then the matrix 

(3-6) [{C\i-t - C-„) ■ ■ ■ {cx,-^ - C^^-J-i)]q<^^^.<^ 

has a block-triangular form with index sets {0, . . . , /c — 1}, {/c}, and {/c + 1, . . . , n}. 
The middle diagonal block is zero and so is the determinant of (3.6). This proves 
(3.5). 

Now let us compute s^(/(i; O). By the same argument the matrix (3.6) with X = fx 
is triangular, therefore 

S/il/^,"J- -pr (c .-c ■) 

= 11 (c^i-i -c^;-j-i) 7 

where one goes from the first line to the second by a standard combinatorial argu- 
ment (used in e.g. proof of the hook- length formula). This proves (3.4). 

There is a convenient tableaux sum formula (see [M2] and [GG]) for the poly- 
nomial (3.3). Given a partition /i, a reverse tableaux T on // is, by definition, a 
function on the diagram of fu. 

^ 3 {i,j) I — > {0,...,n}, 
that decreases along the columns 

T{iJ)>T{i'j), ^>^^ 
and does not increase along the rows 

T(i,j)>T(z,/), j<j'. 
With this notation the tableaux sum formula is 
(3.7) s^(a;;0) = ^ JJ (xT(i,j) - Cj-i-T(i,j)) , 

where the summation ranges over all reverse tableaux T on /i. 
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3.4 Interpolation Macdonald polynomials. 

These polynomials depend on 5 parameters a, 6, c, q, t from^ the field k such that 
q,t y^ and the corresponding grid is given by 

(3.8) [ij]=a + bqH' + ^ek. 
The non-degeneracy of such a grid is assured by inequalities 

q^^t^^ A;>0, />0, 
bqH^^c, /OO, />0, 

see more precise conditions in Appendix A. The simultaneous shift & scale trans- 
formations 

Xi 1-^ CiXi + C2 , i = 0, . . . , n, Ci G k \ 0, C2 G k , 

reduce the number of non-trivial parameters to following 3 

q,t,c/b, 

where the last one can also assume the value cxd. 

In full generality the corresponding interpolation polynomials P^{x; O) were con- 
sidered in [Ok4] . Important particular cases and degenerations of them were con- 
sidered earlier by F. Knop, G. Olshanski, S. Sahi, and the author in a long series 
of papers, see References. 

The following tableaux sum formula for the polynomials P^{x;Q) follows by a 
change of variables from the statement of Theorem 5.2 in [Ok4] 

(3.9) P^{x-9) = Y,^T{q,t-^) n {^n^,J)-[^ + T{^,j),j]) , 

where the summation ranges over all reverse tableaux T on n and 'i(jT{q,t~^) is a 
certain two-parametric weight of a tableaux T introduced by I. Macdonald in [Ml], 
section VI. 7. This weight is a product of factors of the form 

(1-q¥), A;,-/gZ>o, 

and it appears in the tableaux sum formula for the ordinary Macdonald polynomials 

(3.10) P^ix;q,t-^)=J2MQ.t-^) H ^m^) " 



^Actually, the parameters a, b, c, q, t may lie in some extension of the field k; then (3.8) implies 
that q and t satisfy quadratic equations with coefficients in k 
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In (3.10), the LHS denotes the ordinary Macdonald polynomial with parameters 
q and t~^ (we use Macdonald's notation for it; it is not to be confused with our 
interpolation polynomials). Note that, in particular, (3.10) is the highest degree 
term of (3.9). 

The 5-parametric grid (3.8) admits many degenerations; they are all listed in 
Appendix A. As in the two previous examples, the polynomials Pfj.{x; Q) also enjoy 
the property 

P^(A;O)=0, unless /U C A , 

which can be seen from (3.9) (see Lemma 4.1 below) but is actually used in the 
proof of (3.9). 

The polynomials (3.9) seem to be very distinguished special polynomials. They 
have a wealth of applications in various fields of mathematics, see e.g. [Okl,004- 
5,KOO,Ok5]. 

As to the practical interpolation, in the particular case when 

6c = 

there exists a nice efficient algorithm for Newton interpolation, see [Ok3]. It really 
speeds up for the following (Jack) degeneration of the grid (3.8) 

[ij] = a + (3j + (3'i . 

A remarkable feature of that algorithm is that the coefficients of the Newton in- 
terpolation expansion of any symmetric polynomial can be found without actually 
computing these (rather complicated) Newton interpolation polynomials. That 
algorithm can be also used for expansion in ordinary Macdonald and Jack polyno- 
mials. 

Note that in two particular cases 

t = 1 or t = 1/q 

the number (3.8) depends only on j or j — i respectively. Hence for these values 
of t the interpolation Macdonald polynomials become particular cases of factorial 
monomial symmetric functions or factorial Schur functions. 

4. Statement of the characterization theorem 

The exactly solvable examples 2-4 of the previous section have at least two 
following common features. The interpolation polynomials admit a tableaux sum 
formula of the form^ 

(4.1) P^(x;0) = J]weight(T) J] ^XTi^,J) - [i + T{iJ),j]) , 



Ht is easy to see that the formula (3.2) can be written in the form (4.1) with weight(T) £ {0, 1}. 
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where the summation ranges over aU reverse tableaux T on ^ and 

weight (T) G k 

is some weight of a tableaux T. The interpolation polynomials also enjoy the 
following extra vanishing property 

(4.2) P^(A; O) = , unless fi C X . 

In fact, (4.2) is a consequence of (4.1) as the following argument (borrowed from 
[Okl], Section 3.8) shows 

Lemma 4.1. If T is a reverse tableau on n and X is a partition then 

H {[T{iJ),XTi,J~[i + T{iJ)J])=0, 
(*J)eM 
unless fx C X. 

Proof. Suppose that 

[T(i,j),AT(.,,)]-[i + T(z,i),i]^0, V(z,j)G/". 
In particular, for i = we obtain 

(4-3) At(o,o) 7^ , At(o,i) 7^ 1 7 • • • 7 -^T(oj) ¥" J , ■ ■ ■ ■ 

On the other hand, since T is a reverse tableau we have 

(4-4) Xt{o,o) < -^T(o,i) < • ■ ■ ■ 

The inequalities (4.3) and (4.4) imply that 

At(o,j) > i , i = 0, 1, . . . . 

Again, since T is a reverse tableaux we have 

(4.5) T(0,j)<T(l,j)<---<T(^;-l,j) 

and also 

(4-6) j < At(o,j) < At(i,j) < ■ ■ ■ < ^T{^,r-l,J) ■ 

By the definition (2.10) the inequalities (4.5) and (4.6) yield that 

a;- >;«;•, i = o,i,..., 

which is equivalent to /U C A . D 

Remark 4-2. By employing the same argument as used in the proof of Theorem 5.1 
in [Ok4] one can show a priori that, conversely, the extra vanishing (4.2) implies 
existence of a tableaux sum formula (4.1). We shall obtain this implication as a 
corollary of our main theorem. 

Remark 4-3. It is also clear that the extra vanishing (4.2) follows immediately from 
any analog of the g-integral representation established in [Ok4] for the interpolation 
Macdonald polynomials in the case a, 6, c, g, t G C and \q\ < 1. 
The above discussion justifies the following 
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Definition 4.4. We shall call a non-degenerate grid O perfect if the polynomials 
Pfj,{x; O) enjoy the extra vanishing property (4.2). 

Our main result is the following 

Main Theorem. The following is the list of all perfect grids Q: 

El. [ij] = Cj, where cq, ci, . . . are pairwise distinct elements ofk; the correspond- 
ing interpolation polynomials are the factorial monomial symmetric function (see 
Section 3.2). 

E2- [ij] = Cj-i, where . . . , c_i, cq, ci, . . . are pairwise distinct elements ofk; the 
corresponding interpolation polynomials are the factorial Schur function (see Sec- 
tion 3.3). 

I. [ij] = a + hqH'^ + cq~H~\ where a,h,c,q,t are elements of a certain extension 
of k; the corresponding interpolation polynomials are the interpolation Macdonald 
polynomials (see Section 3.4). 

II-IV. The grid O and the interpolation polynomials are one of the degeneration of 
the previous case (see Appendix A). 

Since the two first cases are much simpler than the remaining ones we refer to 
them as to the 1st and 2nd elementary cases and number them by Ei and E2. 

It follows from the above theorem together with Lemma 4.1 and Remark 4.3 that 
if there exist any other exactly solvable cases of the general interpolation problem 
described in Section 2, then the formulas for the corresponding interpolation poly- 
nomials should have some entirely new structure. 

The proof of the theorem will be given in Sections 5-7. 

4. Reductions of the proof 
First from (2.3) and (2.9) one immediately derives the following 
Proposition 5.1. If a grid O 

(5.1) O : {0, ...,?i} X Z>o — >k 

is perfect then so are all grids Q,^ , /c = 0, 1, . . . , and O^, to = 0, . . . , n. 

Introduce one more operation on grids. Given a grid (5.1) we can define a grid 

;0 : {0, . . . ,n — /} — >k, / = 0, . . . , n , 

by setting 

zO(z,i) :=0(z + /,i). 

From Proposition 2.3 it is clear that the operation 

O I — > iQ 

preserves non- degeneracy. We now plan to show that it preserves perfectness as 
well. First, we establish the following 
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Proposition 5.2. Let Q be a perfect grid and let P^{x;Q) be the corresponding 
Newton interpolation polynomials. We have 

(5.2) P^(a:;0)=x'^ + ..., 

where x^ = Xq° ■ ■ ■ x^" and dots stand for lower monomials in lexicographic order. 
Proof. First show that 

(5.3) P^(x;0) = c^x^ + ... , c^Gk, 

for certain constants c^. Induct on \fx\, the case |/u| = being clear. Let m 

m:=deg^^^Pf,{x;n) 

be the degree of P^(x; O) as ofa polynomial in xq. If m < /Uq then (5.3) is established 
with c^ = 0. 

Suppose therefore that m> \no\ and consider the leading coefficient 

g{xi, ...,Xn) := [x^] P^ix] O) G k[xi, . . . , x. 



lS(n) 
L, . . . , ^n\ 



of P^{x] O) as of a polynomial in xq. By our hypothesis g is a, non-zero polynomial 
of degree 

degfif < \^\ - ^0 = Hi -\ h //n • 

We claim that 

(7([l,Ao],...,[nA^_i])=0 

for all partitions A such that A^ < fXi+i for some i = 0, . . . , n — 1. Indeed, by the 
extra vanishing condition (4.2) the polynomial 

Pf, {xo, [1, Ao], . . . , [nA„_i]; O) 

has in this case infinitely many zeros 

xo = [0,Ao],[0,Ao + l],[0,Ao + 2],... 

and hence vanishes identically. Since the grid lO in non-degenerate we conclude 
that 

g = const P(^^,...,^„)(xi,...,a:n;iO) , 

which by inductive assumption establishes (5.3). 

Now recall that for any non-degenerate grid the polynomials P^{x;fl) form a 
linear basis in the k-linear space of all symmetric polynomials. This immediately 
implies that 
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Let us renormalize the polynomials P^{x; O). Introduce new polynomials 

P^(a;;0):=— P^(a;;0) 

for which we have 

P^(x;0)=x'^ + .... 

It is clear that these polynomials also satisfy the equalities (2.3) and (2.9). There- 
fore, by Proposition 2.9 we conclude that 

P^(x;0)=P^(x;0), 

which concludes the proof of the proposition. D 

In fact, the above argument establishes more than just (5.2). We have also 
proved the two following facts: 

Proposition 5.3. We have 

P^{x; O) = a:({°P(^^,...,^„)(a;i, . . . , x^; lO) + . . . , 

where dots stand for terms of lower degree in xq. 

Proposition 5.4. // a grid O is perfect then so are the grids lO, 2^^, • • • • 
Let us combine the 3 operations 

as follows. Given a grid (5.1) we define for all < / < m < n and all /c > a new 
grid 

zO^:{0,...,m-/}xZ>o-^k 

by the formula 

zO^(z,i):=0(z + /,i + /c). 

Then by Propositions 5.1 and 5.4 we have 

Proposition 5.5. // a grid O is perfect then so are all grids zOj^. 

From now on we assume that O is a certain given perfect grid and our goal is to 
show that O is one of the grids listed in the Appendix A. 

Introduce the following notation. By the symbol e we shall denote some non-zero 
element of k 

e Gk\0. 

The purpose of this notation is to denote irrelevant overall factors in our equations. 
An example of an element e is any product of factors of the form 

Now we prove the following 
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Proposition 5.6. Let F he the following rational function 

(^ .. T-., . yl ~ y4yi - y4y5 + y2y3 + yiy5 - Vs 

(5.4) F(yi,y2,y3,2/4,2/5) := • 

?/2-?/3 

Then we have 

(5.5) [i + l,j + 2] = F{[i,j],[i + l,j],[i,j + l],[i + l,j + l],[i,j + 2]), 
for all j and all i < n — 1 . We also have 

(5.6) [i + 2, J + 1] = F{[i,j], [i,j + 1], [z + 1, j], [i + l,j + 1], [i + 2,j]) , 

for all j and all i < n — 2. 

Proof. Prove (5.5). By virtue of Proposition 5.5 it suffices to establish (5.5) for 

n = l, i = j = . 

Consider the condition 

(5.7) P(3,o)(('?2);n)=0. 
A direct computation shows that 

(5.8) P(3,o)((2r2);fi) = 

= £(HN - [00][ll] - [of + [0l][l0] + [0l][l2] - [02] [ll] - [10][12] + [if) , 

where 

_ ([12]-[11])([12]-[10]) 
[01] - [10] 

Since 

[oi] - [lo] ^ 

we deduce from (5.7) that 

(5.9) [i2]=F([oo],[io],[oi],[ii],[o2]). 

This proves (5.5). 

Prove (5.6). By virtue of Proposition 5.5 it suffices to estabfish it for 

n = 2, i = j = . 

Consider the condition 

P(2,o,o)((iTM);fi) = o. 
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A direct computation shows that 

(5.10) P(2,o,o)((iTm);0) = 

= £([10][21] - [01] [21] + [if - [ll][00] + [0l][l0] - [10]2 + [00][20] - [ll][2o]) 

where 

^ [21] - [20] 



[10] - [01] • 
Therefore, 

(5.11) [2l]=F([00],[0l],[l0],[ll],[20]). 

This concludes the proof. D 

It is clear from the formula (5.4) and non-degeneracy that the equalities (5.5) 
and (5.6) can be reversed as follows 

Corollary 5.7. We have 

(5.12) [ij] = F{[i + l,j + 2], [i,j + 2], [z + 1, J + 1], [i,j + 1], [z + 1, j]) , 
for all j and i < n — 1, and also 

(5.13) [ij]=F{[i + 2,j + l],[t + 2,j],[i + l,j + l],[t + l,j],[i,j + l]), 

for all j and i < n — 2. 

The above equalities immediately result in the following two propositions. 
Proposition 5.8. Suppose that for some i and j we have 

[i,j] = [i + l,j + l]. 

Then Q is a grid of type E2 that is, [z, j] depends on j — i only. 

Proof. Follows from (5.5), (5.6), (5.12), (5.13), and the following identity 

(5.14) F{z,u,v,z,w) = V . D 

Proposition 5.9. Suppose that n = 1 and for some j we have 

[0,j] = [l,j], [0,j + l] = [l,j + l]. 
Then Q is a grid of type Ei that is, for all k > we have 

[0, k] = [1, k] . 

Proof. Follows from (5.5), (5.12), and the following identity 

(5.15) F{u,UjVjV,w) — w . n 
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6. Proof of the Theorem in the two variables {n = 1) case. 
Consider the condition 

P(4,o)((3^);0) = o 

and make the substitution (5.11). After that substitution one computes: 

(6.1) P(4,o)((1U); fi) = e (H - [ii]) Fi (H, [lo], [oi], [ii], [02], [03]) , 
where e is the following invertible element 

_ ([03]-[02])([l2]-[ll])([l2]-[l0]) 
([02] - [10])([10] - [0l])2 

and Fi denotes the following polynomial 

Fi{yi, ?/2, 1/3, ?/4, 1/5, ye) ■■= ?/6G'i(?/2, 2/3, 2/4) - Go{yi, 2/2, 2/3, 2/4, 2/5) , 
where 

Go(2/i, 2/2, 2/3, 2/4, 2/5) := 2/12/42/2 - 2/12/42/5 - 2/12/22/5 + 2/12/5 + 2/32/4 + 2/32/2 

- 2/32/5 - 2/32/4 - 2/32/42/2 - 2/32/2 + 2/32/5 + 2/42/5 + 2/42/22/5 - 2/42/5 + 2/22/5 - 2/22/5 

and 

^1(2/2, 2/3, 2/4) := (2/3 - 2/2) (2/3 - 2/4) • 

Since the product (6.1) vanishes we have the following alternatives. 

Case 1. We have [00] — [11] = 0. Then by Corollary we find ourselves in the E2 
case. 

Case 2. We have 

(6.2) [o3]Gi([io], [01], [11]) - Go([oo], [10], [01], [11], [02]) =0. 

Then, since the above equation is linear in [03], we have again two possibilities. 
Case 2.1. We have 

(6.3) Go ([00], [10], [01], [11], [02]) = Gi([io], [01], [11]) =0. 
Then the only solution of 

Gi([io],[oi],[ii])=([oi]-[io])([oi]-[ii]) = 
compatible with the non-degeneracy is 

(6.4) [01] = [11] . 
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Substituting it into the first equation we obtain 

([02]-[l0])([02]-[0l])(H-[l0])=0. 

From non-degeneracy we obtain 

[oo] = [lo] . 

Thus, by Proposition 5.9 we are in the Ei case. 

Case 2.2. We assume that the grid O is not of the elementary types Ei or E2. 
Then we can uniquely determine [03] from the values of [00], [10], [01], [11], and [02] 
by the formula 

G'o([oo],[io],[oi],[ii],[o2]) 



[03] = 



Gr{[io],[oi],[ii]) 



Then using (5.5) we can also uniquely determine [13]. We claim that we can deter- 
mine [04] as well. To that end we repeat the entire argument for the grid O^ which 
is also perfect by Proposition 5.5. We claim that for the grid fl^ the only possible 
case is again the Case 2.2. Indeed, the grid O^ cannot be of types Ei or E2 because 
otherwise by Propositions 5.8, 5.9 the grid O itself is of type Ei or E2. Therefore, 
we can uniquely determine [04] and then [14] and then all the rest 

[ij] , z = 0, 1 , j>5. 

in the same manner. 

Now, this implies that our grid O is one of the types I-IV. Namely, we shall check 
below in Proposition 6.1 that for any perfect grid O which is not of the elementary 
types El or E2 one can find a grid O of one of types I-IV such that 

[00] = [00] , [10] = [10] , [01] = [01] , [11] = [11] , [02] = [02] . 

But then, since those five values uniquely determine the rest, we conclude that 
= 0. 

Therefore, to finish the proof of the theorem in the n = 1 case it suffices to 
establish the following 

Proposition 6.1. Denote by E the following subset o/Z>o x Z>o 

E = {(0, 0), (1, 0), (0, 1), (1, 1), (0, 2), (0, 3)} C Z>o x Z>o . 

Then any perfect grid O belongs to one of the following classes (compare with clas- 
sification in Appendix A): 

El. [Oj] = [lj],j = 0,1,2,.... 

E2. [0,j] = [l,j + l],j = 0,l,2,.... 
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I. There exist a, b, c,q,t Ek, q^t ^ 0, ±1 such that 

[ij]=a + 6Q¥ + ^, V(i,j)GE. 

The elements a, 6, c, g, t G k are determined uniquely up to the following symmetry 

q^^q~^, t h^ t~^ , h^^c, cy-^h. 

II. There exist unique a, /?, /?', 7 G k such that 

[ij] =a + (3j+ (3't + -f{(3j + (3' if , V(z, j) G S . 

Ilia. There exist unique a, a', 13, 13' G k snc/i that 

[ij] = « + (-1)^'(«' + (3j + (3'i) , V(z, j) G S . 

Illb. There exist a, a', /?, /?' G k snc/i t/iat 

/n ^/izs case on/y the numbers a + a' , a — a' + (3' , and (3 are determined by the 
numbers [ij], {i,j) G E. 

IIIc. There exist unique a, a', /?, /?' G k snc/i t/iai 

[ij] =a + {-lY+^{a' + /3j + /3'z) , V(i,i) G E . 

IV. There exist unique a, /?, /?', g G k snc/i that 

[Qj]=a + I3q\ j = 0,...,3, 
[li]=a + /3V^ i = 0,l. 

The proof will be based on the following lemma which can be established by 
direct inspection 

Lemma 6.2. For any 4-tuple 

{wo,Wi,W2,W3) G k"^ 

satisfying 

Wi ^ W2 

we have the 3 following mutually exclusive possibilities: 
(i) There exist a, 6, c, g G k, q j^ 0, ±1 such that 

Wj = a + bq^ + — , j = 0, . . . , 3 . 
qj 
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The elements a,b,c,q Ek are determined uniquely up to the following sym- 
metry 

q^^q~^, bi-^Cj Ci-^b. 

In this case 

Wo — 3wi + 3w2 — ws y^ and wq + wi — W2 — ws j^ . 
(ii) t^o — 3tyi + 3w2 — ws = and there exist unique a, /3, 7 G k such that 

Wj = a + Pj + 7/?^/ , j = 0, . . . , 3 . 
(iii) wq + wi — W2 — ws = and there exist unique a, a', /3 G k such that 
[ij]=a+i-iyia' + Pj), j = 0,...,3. 

The last case is a subcase of second one z/char k = 2 . 

Proof of Proposition 6.1. Assume that we are not in the elementary cases Ei or E2 
and apply Lemma 6.2 to the following 4-tuple 

[00], [01], [02], [03] G k. 

Recall that by non-degeneracy we have [01] 7^ [02]. Therefore we have 3 mutually 
exclusive possibilities which we shall call cases (i), (ii), and (iii) respectively. 

(i) In this case we have 

(6.5) [0i]=a + 6Q^ + 4' J = 0,...,3. 

Assume for a moment that be 7^ 0. Choose to ^ k so that 



(6.6) [10] = a + 6to + 



c 



Then the other root of this equation equals - — . Since we exclude the elementary 

bto 
cases the following equation is satisfied and is not identically zero 

(6.7) [o3]Gi([io], [01], [11]) - Go([oo], [10], [01], [11], [02]) =0. 

Substituting (6.5) and (6.6) into (6.7) we obtain 



q^b-c) (a + bqto + — - [ 
V Qto 



nll(a+^ + ^-Nl=0. 
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Since [oi] 7^ [02] we have bq'^ 7^ c. Hence either 

11 =a + bqto + -— , 

qto 

in which case we are done by setting t = to, or 

r 1 bto cq 

11 =a^ + — , 

q to 
1-1 c 

m which case we set t = — . 

bto 
Now consider the case be = 0. We can assume that c = 0. Set 

[10] - g 
b • 

Then the equation (6.7) imphes that either 

[11] = a + bqt , 
which means we are in case I, or 

r 1 bt 

11 = a-\ , 

Q 

which brings us in the case IV. 

(ii) In this case we have 

(6.8) [Oj]=a + /3j+7/3V, j = 0,...,3. 

By non-degeneracy we have P ^ 0. Again, assume for a moment that 7 7^ 0. Choose 
(3q Ek so that 

(6.9) [io]=a + /3^+7/3^'. 

Then the other root of this equation equals —(3q — -. Solving the equation (6.7) 
for [11] we obtain that either 

[ii]=a + /? + /3^ + 7(/3 + /3^)2, 

in which case we set (3' = (3q, or 

[ii]=a- 13 + (3'o +-f {-(3 + (3'o)\ 

in which case we set /?' = —(3q — - to obtain 

[11] = a + /3 + /3' + ^{/3 + /3'f. 

Now if 7 = we set (3' = [10] — a and the two possibilities for [11] bring us in the 
cases II and Illb respectively. 

Similarly, the consideration of the (iii) case leads to cases Ilia and IIIc. This 
concludes the proof. D 
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7. Proof of the theorem for n > 1 
Suppose O is not of type E2; then by Proposition 5.8 we can assume 

(7.1) [00] ^ [11] . 
Consider the condition 

(7.2) P(3,o,o)((?M);0) = 0. 
Using the substitution (5.11) one computes 

(7.3) P(3,o,o) ((27571);^) = 

- £([11] - H) ([20] G2([oo], N, [01], [11], [02]) - G3([oo], [10], [01], [11], [02])) 

where 

([02]-[0l])([2l]-[20])([ll]-[20]) 



e 



([20]-[0l])([l0]-[02])([l0]-[0l])2 

and G2 and G3 are the following polynomials 

(7.4) G'2(?/i,y2,y3, 1/4,1/5) ■=y4y5 -1/41/3 + 2/21/1 -1/21/3 -yiy^ + vh 

Gsivi, ?/2, 1/3, ?/4, 1/5) := yl - yly2 - yjyi - ?/|i/3 - y^yl + y^y2yi + y4y2y5+ 

22/42/22/3 + 2/42/12/3 - 2/42/3 + 2/1 - 2/12/5 - 2/12/3 + 2/22/52/3 - 2/22/3 - 2/i2/32/5 + 2/| • 

Recall that we assume that [00] 7^ [11]. Therefore, from (7.2) we have a linear 
equation in [20] and hence two possible cases: 

Case 1. We have 

^2 = ^3 = . 

Using (7.1) we can express [02] from G2 = 

^ [01] ^-[10] [01] -[11] [01] + [10] [00] 
[00] - [11] 

Substituting this expression into the equation G3 = we obtain 

(7.6) £([11] - [oi])([ii] - [00] - [01] + [io])([ii] - [00] + [01] - [10]) = 0, 

where e = [11] — [10]. Observe that it is impossible to have [11] = [01] because then 
from (7.5) we conclude [02] = [10], which contradicts non-degeneracy. Similarly, it 
is impossible to have 

[11] = [00] + [01] - [10] 
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because then from (7.5) we conclude 

[02] = [00] 

again in contradiction with non- degeneracy. Thus, the only possible solution of 

(7.6) and (7.5) is 

[11] = [00] - [01] + [10] , 
[02] = 2 [01] - [00] . 

That is, the grid O2 is of the type Illb. Recall that the parameters a, a', /3, 13' are 
not uniquely determined by the grid O2 (see Proposition 6.1). We can find such 
values of a, a', /?, /?' G k that 

[ij] = « + (-l)^(a' + (3j + (3'i) , i+j<2. 

for all « + J < 2. Then by (5.6) we conclude that 

(7.7) [ij] = a + {-iy{a' + I3j + I3'i), Vi < 2 V j > . 

Note that now the parameters a, a', /3, [3' are uniquely determined. If n = 2 this 
finishes the consideration of the Case 1. 

If n > 3 then we have to look at the equation 

(7.8) P(2,i,o,o)((lXiri))=0. 
By (7.7) one concludes from the identity 

F(w, u + z^v^v — z^w) = w + z 
and (5.6) that 

(7.9) [31] = [30]-/?. 

From (7.7) and (7.9) one computes that (7.8) is equivalent to 

{a-a' - 3(3' - [30]) ([21] - [30]) = 
By non-degeneracy we conclude that 

[30] = a — a' — 3(3' 
and then by (5.6) and (7.7) it follows that 

[3j] =a-a' -3(3' - (3j . 
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Thus, the grid O is of the type Illb. Repeating the above argument for the grids 
lO, 2^^, ... we conclude by Proposition 5.5 that the entire grid O is of the type Illb. 

Case 2. Now we assume that Q is not of the Illb type and then we can determine 
[20] from (7.3) 

(7.10) [20] = g3([oo],[io],[oi],[n],N) _ 

G'2(HJl0]j0l]Jll]j02]) 

After that using (5.6) we can also determine [21], [22], .... It is clear that if the grid 
Oi is of the type Ei,I,II,IIIa,IIIc then the grid Q2 is of the same type with the same 
parameters. If the grid Oi is of type IV then 

^2 = , Gs ^ , 

that is, the equation (7.2) has no solutions. For n = 2 this concludes the proof. 

If n > 2 then we have to repeat the entire argument for the grids lO, 2^, ■ ■ ■ ■ In 
order to be able to do so we have to verify that once the grid O2 is not of type Illb 
then neither is the grid 1O3. But this follows from the just established classification 
of the perfect grids for n = 2. Namely, it follows that if O is a perfect grid then 

(1O2 is of type Illb) =^ (O2 is of type Illb) . 

This concludes the proof of the theorem. 

Appendix A. Table of perfect grids. 
El. The first elementary case 

where 70,71,... are arbitrary pairwise distinct. The corresponding interpolation 
polynomials are factorial monomial symmetric functions. 

E2. The second elementary case 

M = ij-i , 

where . . . , 7_i, 70, 71, . . . are arbitrary pairwise distinct. The corresponding inter- 
polation polynomials are factorial Schur functions. 

I. The generic case of Macdonald interpolation polynomials 

\ij] = a + hqh'^ + -^ , 

where g, t 7^ 0. The non-degeneracy conditions are 

q'' ^t\ V/c, / /c> 0, < / < n , 
bqH^^c, yk,l k>0,0<l<2n. 
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If t = 1 or t = 1/g we hit the two elementary cases above. 
II. The case 

can be obtained from I by setting 

7 1 7 1 7 

q = l+/3h, t = l+/3'h 
and letting h ^ 0. The non-degeneracy conditions are 

k(3^l(3' , \/k,l k>0,0<l<n, 
-f{k(3 + l(3')^-l, ykj /c>0,0</<2n, 

The cases II and Ill(abc) are impossible if char k > . 
Ill(abc). The case 

can be obtained from I by setting 



a' 1 _ a' _ 1 

q^e{l+/3h), t = e{l+/3'h) 



a = a , 0=^^ + 77;-, c = — — 77;-, 



and letting /i ^ 0. Here e, e' = ±1; namely, we have three subcases: (a) (e, e') = 
(-1, 1), (b) (e, e') = (1, -1), (c) (e, e') = (-1, -1). 

IV. This case exists only for n — 2 when we can have 

[Oj] = a + (3q\ [lj]=a + (3'q-^. 

It is obtained from I by setting 

a = a , h = (3 , c = 13' t 

and then letting t -^ 0. The non-degeneracy in this case is equivalent to 13,(3' ^ Q 
and q^ ^ 1,13'/ 13 for aU /c> 0. 
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Appendix B. An example of a universal interpolation polynomial. 

Consider the case n = 1 and 

/" = (3,0). 

Then the universal polynomial P'^{xq,xi) lies in 

1 1 



P^(xo,Xi)GZ 



Uqo,Uqi,Uo2,Uio,Uii, 



[xq.Xi]. 



U02-U10 Uoi-Uio_ 

In the basis of the monomial symmetric functions m\ it is given by the formula 

™/ N , C2i(w)m2i + 020(^)^20 + cii('u) mil + cio(w)mio + coo(w) 
-r„ [xo, xi) — rriso H — , 

[Uo2-Uio)[Uoi -Uio) 

where C21, . . . , cqo are the following polynomials 

C2I =Uo2Uoo - U02U10 + U02U01 - U02U11 - UooUio + UooUoi - UooUii + uIq 

- uiouoi + uiouii - uoiuii + ul^ , 

C20 ='"02^10 - UI2U01 - U02U00U01 + U02U10U01 + U02U10U11 - Uo2uli + UooUioUii 

- -ufo - uIqUii + uiouli + WloWoi^ll - WlO^ll 5 

Cll = - {uq2 + Uio + Uoi + Uii){uo2Uoo - U02U10 + UQ2U01 - UQ2U11 - UoqUiq 

+ UoqUqi - UooUii + uIq - UiqUqi + UioUii - Woi^ll + ^ll) , 
ClO =W02^00W01 - f^02«10 - UI2U10U11 + ^02^01 + f^02WooWoi + Uo2ulo - Uo2uloUoi 

- U02U10U01U11 - UoquIqUii - UooUiouli + uIqUqi + uIqUii - Wio^oi 
+ ^10^11 - WlO^oi^ll + ^10^11 , 

Coo = - U02^00uli + ^02^10^01 + ^02^10^11 " ^02^10^01 - Uo2uIqUoi - Uo2uIqUii 

+ uo2ulouli + uo2uIqUoiUii + Woo^io^ii " ^lo'^oiWii + ^lo'^oi^ii - '"10^11 • 
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